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Critical current density was studied in single crystals of Ba(Fei_ :E Co :c )2As2 for the values of 
x spanning the entire doping phase diagram. A noticeable enhancement was found for slightly 
underdoped crystals with the peak at x = 0.058. Using a combination of polarized-light imaging, 
x-ray diffraction and magnetic measurements we associate this behavior with the intrinsic pinning 
on structural domains in the orthorhombic phase. Domain walls extend throughout the sample 
thickness in the direction of vortices and act as extended pinning centers. With the increasing 
x domain structure becomes more intertwined and fine due to a decrease of the orthorhombic 
distortion. This results in the energy landscape with maze-like spatial modulations favorable for 
pinning. This finding shows that iron-based pnictide superconductors, characterized by high values 
of the transition temperature, high upper critical fields, and low anisotropy may intrinsically have 
relatively high critical current densities. 

PACS numbers: 74.70.Dd,72.15.-v,74.25.Jb 



INTRODUCTION 



In type II superconductors magnetic field penetrates 
the bulk in the form of Abrikosov vortices. In the 
presence of electric current these vortices experience the 
Lorenz force, and, if not hindered, their motion leads to 
a non-zero resistance. Vortices can retain their positions 
due to pinning caused by local variations of the supercon- 
ducting properties that lead to a position-dependent vor- 
tex energy. Therefore, electric current flows without dis- 
sipation, as long as its density, j, stays below some tem- 
perature and field dependent critical value, called criti- 
cal current density, j c The pinning strength is low 
in homogeneous high-quality single crystals, unless they 
have defects introduced either artificially or naturally [2j . 
There are different types of pinning centers ranging in 
size and dimensionality. Since vortices are linear ob- 
jects threading the entire sample, they are more strongly 
pinned by extended defects. In highly anisotropic lay- 
ered superconductors pinning is maximum when vortices 
lay parallel to the layers. However, the measured criti- 
cal current density is determined by a much weaker pin- 
ning of vortices that are perpendicular to the supercon- 
ducting layers The most efficient way to stop such 
vortices from bending and moving is to pin them by 
defects that extend perpendicular to the layers. Two- 
dimensional (planar) pinning centers can appear natu- 
rally, for example, as twin boundaries in YBa 2 Cu307_ y 
(YBCO) or in the orthorhombic/antiferromagnetic 
(AFM) phase in ErNi 2 B 2 C and HoNi 2 B 2 C @, S|. Lin- 
ear pinning centers can be artificially created, for ex- 
ample, by heavy ion irradiation that produces colum- 
nar defects (see, e.g., Section IX in [2]). Whereas such 
defects lead to certain enhancement of pinning in high 



critical temperature cuprate superconductors (high-T c 
cuprates), their extreme anisotropy has been proven to 
be detrimental for technological use. The discovery of 
iron-based pnictide superconductors @, 0] with relatively 
high transition temperatures, T c , has inspired new hopes. 
Not only do these materials have very high upper criti- 
cal fields comparable to cuprates [1] , they have relatively 
low anisotropy 0, E3, EH 0, EH , which is beneficial for 
applications. Nevertheless, the vortex properties of the 
iron pnictide superconductors were found to be similar 
to the cuprates and can be understood within the weak 

collective pinning approach 14|, [l 1E1,E3,E3,E1- Not 

surprisingly, heavy ion irradiation of pnictides also leads 
to the enhancement of the critical current density [2(| • 



In this work, we show that by fine tuning the composi- 
tion in Ba(Fei_ a; Co a ;) 2 As 2 ("Co-Bal22") one can achieve 
a significant increase of pinning due to twin boundaries 
in the AFM/orthorhombic phase. Previously we directly 
imaged such boundaries in several AEFe 2 As 2 (AE=Ca, 
Sr, Ba) parent compounds 211 ] . The effect of domain 
boundaries on the direction of vortex motion in Co- 
doped compositions was shown with SQUID microscopy 
[22J. Here we demonstrate that these domains co-exist 
with superconductivity almost up to optimal doping and 
play the role of efficient extended pinning centers on the 
slightly underdoped side. Unlike YBCO, where the den- 
sity of twin boundaries varies with sample and, in the 
cleanest samples, can be absent, in the pnictide single 
crystals a maze of fine robust domains always appears 
up to the optimal doping level and leads to a substantial, 
intrinsic pinning. 
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EXPERIMENTAL 

Single crystals of Ba(Fei_ a; Co a ;)2As2 were grown from 
FeAs flux from a starting load of metallic Ba and FeAs 
and Co As, as described in detail elsewhere Q. Crystals 
were thick platelets with sizes as big as 12x8x1 mm 3 
and large faces corresponding to the tetragonal (001) 
plane. Samples for magnetization and optical microscopy 
polarized-light imaging were cleaved with a razor blade 
into rectangular shaped platelets of typical dimensions, 
(2 - 5) x (2 - 5) x (0.1 - 0.5) mm 3 . 

To characterize the lattice chaqnge associated with 
the tetragonal-to-orthorhombic structural phase transi- 
tion and to check the crystal perfection, high-energy 
x-ray diffraction measurements (E = 99.54 keV; A = 
0.01246 nm), using an area detector positioned 1597 mm 
behind the sample, were performed on the high energy 
station (6ID-D) in the MUCAT Sector at the Advanced 
Photon Source. At this high energy, x-rays probe the 
bulk of the sample rather than just the near-surface re- 
gion and, by rocking the crystal about both the hori- 
zontal and vertical axes an extended range of a chosen 
reciprocal plane can be recorded [23j]. For these measure- 
ments, the samples were mounted on the cold-finger of a 
closed-cycle refrigerator using Kapton windows to avoid 
extraneous reflections associated with beryllium or the 
aluminum housing. For each data set, the sample was 
rotated over a range of horizontal angles, fi, of ±3.2 deg 
for each value of the vertical angle, 77, between ±3.2 deg 
with a step size of 0.4 deg. The splitting of the tetrag- 
onal (220) T diffraction peak into the (400) o /(040) o or- 
thorhombic diffraction peaks signals the tetragonal to or- 
thorhombic distortion. The measured splitting was used 
to determine the dimensionless orthorhombicity param- 
eter 5 = (a - b a ) I (a a + b ) : where a and b are the 
in-plane lattice constants of the orthorhombic unit cell. 

Polarized-light imaging of domain structure was per- 
formed in a flow-type liquid 4 He cryostat with the sample 
in vacuum. The sample was positioned on top of a pol- 
ished copper cold finger and directly observed under lin- 
early polarized light through an analyzer. The analyzer 
was turned, close but not exactly at 90° with respect to 
a polarizer. Due to the anisotropy of the optical proper- 
ties in the orthorhombic phase, there is a rotation of the 
polarization plane upon reflection, always from a— to b— 
axis and the amount of the rotation scales with the or- 
thorhombicity parameter 5. Thus we observe a contrast 
between adjacent domains, because they rotate light in 
opposite directions. 

The magnetization measurements were conducted in a 
Quantum Design MPMS magnetometer. In the experi- 
ment, a platelet sample was fixed in a gelatine capsule 
with a small amount of Apiezon grease and the capsule 
was placed in a clear plastic straw. Care was taken to 
avoid artifact related to the tilt of the samples by re- 



peated removal, re-mounting and re-measuring the same 
sample. At some concentrations several samples of dif- 
ferent aspect ratio were measured. Then the measured 
magnetization was converted into the critical current 
density producing practically identical values indicating 
good quality of the samples and applicability of the Bean 
model. Critical current densities were determined from 
the measured magnetization following the Bean model 
[H , in which for a rectangular slab with the dimensions of 
2c < 2a < 26, j c [A/cm 2 ] = 20M[emu]/ [aV(l ~ a /(3b))} 
where all dimensions are in cm, M is the total measured 
magnetic moment and V is the sample volume. The ob- 
tained values are in a good agreement with those deter- 
mined from the measurements of the magnetic induction 
profiles as well as from direct transport measurements as 
shown previously 15j. 

The critical current density was compared at T/T c = 
0.7 that was high enough to reach the field of the fish- 
tail magnetization maximum, H m , see Fig. El which has 
been used as a characteristic field. Another possibility 
was to evaluate the current in the remanent state, but 
the value of a magnetic moment at H =0 is significantly 
affected by the sample geometry due to demagnetization. 
Evaluation of the current density at H m has avoided such 
problem. More details on the measurements and conver- 
sion procedures can be found elsewhere [Hj]. 



RESULTS 

In the previous study of domain formation in the par- 
ent compounds BaFe2As2, we showed that the domains 
are of 45° type, i.e. their walls are along the (110) direc- 
tion in orthorhombic notation (or parallel to the tetrag- 
onal (100) axis) and that they span along the c-axis 
through a significant fraction of the crystal of 0.2 mm 
thickness. In Fig. Q] we show the evolution of the do- 
main pattern with Co doping, x. As x increases, several 
changes occur. (1) The characteristic separation between 
individual domain walls gets smaller, Fig. [IJd), so we 
could not obtain a digital image that resolves individual 
boundaries already for x > 0.024. However, we could still 
see them by eye up to x < 0.047, with the resolution close 
to the optical limit ~ l/j,m set by the wavelength. On the 
other hand, the domains are large enough so that they 
do not cause X-ray Bragg diffraction spots to broaden 
above the instrumental resolution limit of about 150 nm. 
This gives an estimate of the characteristic width of the 
structural domains between 0.2 and 1 /im for x=0.047. 

(2) As a result, the density of the domain walls increases. 

(3) The same three typical domain patterns are found in 
parent and doped compounds, with vertical, horizontal 
and interwoven (crossing) boundaries. The interwoven 
patterns of domain bundles (not resolving the individ- 
ual boundaries), shown in Fig.[ljc) for x=0.054, become 
much more frequent with increasing x. (4) Optical con- 
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FIG. 1: Polarized-light images showing structural domains 
at 5 K in Ba(Fei_ I Co a; )2As2 crystals for a;=0 (a), 0.0135 (b) 
and 0.054 (c). With increasing x, the domain wall separation 
d becomes smaller (see intensity profiles across the boundaries 
in bottom right panel (d), giving d «4 fim for :r=0 and d ~3 
/im for £=0.135. For x >0.024 we could not resolve individual 
domains since their width is below our digital camera reso- 
lution limit of ~2 fim, however, the domains could be still 
tracked by eye. With further doping, the contrast of individ- 
ual boundaries fades away and they become indistinguishable 
for a;=0.047. However, wide stripes of domains ensembles can 
be seen for x up to 0.054. Eventually the pattern becomes 
optically indistinguishable for x > 0.058 (not shown here). 



trast of domain boundaries fades away with an increase 
of x. 

We could not optically resolve either individual do- 
main boundaries or domain bundles in the compositions 
with x >0.054, despite the observations of anomalies in 
the temperature dependence of resistivity p(T) accom- 
panying structural transition at T s > T c for sample with 
a;=0.058 Q and a notable suppression of the low temper- 
ature phonon part of thermal conductivity at tempera- 
tures below 1 K due to scattering off domain boundaries 
in samples with x up to 0.074 [24J. The difficulty in re- 
solving domains for x close to optimal doping is primar- 
ily due to a suppression of the degree of the orthorhom- 
bic distortion, characterized by dimensionless parameter 
8{x). The doping dependence of 8 is shown in the inset 
in top panel of Fig. [2j The decreasing 8 affects the opti- 
cal contrast, since linearly-polarized light imaging detects 
the anisotropic optical response of the neighboring twins 
[lH and the degree of this anisotropy is proportional to 
8. This decreases the difference in the light intensity be- 
tween the twins. 

As can be seen from top panel of Fig. [5] lines of mag- 
netic, T m (x), and structural, T s (x), transitions coincide 



in the parent compound, x = 0, and then separate with 
the increasing x j^, [26| . These lines can be extrapolated 
to T = to estimate the critical concentrations. Keeping 
in mind that the experimental data are confined to high 
temperatures with T Sl T m > T c and thus the extrapola- 
tion is not precise, we estimate the critical concentration 
for the structural transition, x s «0.070±0.03, and for the 
magnetic transition, x m « 0.065±0.03, see dashed lines 
in the top panel of Fig. O This opens a window of con- 
centrations between x m and x s where only the structural 
distortion is present. The X-ray measurements of the or- 
thorhombic distortion indicate that 8(x), taken at a base 
temperature of 7 K, also decreases and extrapolates to 
zero in the same x range within the experimental error, 
see the inset in Fig. [2l Importantly, since magnetic order 
in 122 pnictides has preferential direction in the ab-plane 
27| and is rigidly correlated with the orthorhombic axes, 
structural domain walls are accompanied by the disrup- 
tion of the magnetic order. In other words, magnetic 
domains follow the structural domains. A similar sit- 
uation of correlated structural/magnetic distortion in a 
magnetically ordered phase is found in the borocarbide 
superconductor ErN^I^C 0, H||. Despite the uncer- 
tainty in the determination of the exact positions of x m 
and x s , they are clearly quite close to the optimal con- 
centration, so that superconductivity coexists with both 
types of order. 

In the bottom panel of Fig. [2] we show the evolution of 
the superconducting critical current density j c .max, taken 
for all doping levels at the same reduced temperature 
T/T c (x)=0.7 at a magnetic field corresponding to the 
maximum in the j c {H) "fishtail" curve, H mi as shown 
by the arrow for x =0.061 in Fig. [3] The magnetization 
curves are shown in the units of current density to allow 
for a direct comparison between samples of different sizes. 
Details of the conversion based on the Bean model [ij for 
finite samples are given elsewhere [Hj]. There is a clear 
asymmetry of j c ,max{x) with respect to the T c (x) dome. 
As a function of doping, the critical current density starts 
to rise with decreasing x below approximately x s , peaks 
at x m or slightly below, and drops off quickly on mov- 
ing towards x = 0. The observed doping dependence of 
the critical current suggests that the best pinning condi- 
tions are realized in the presence of structural domains, 
and, perhaps, with the additional contribution of mag- 
netic pinning. We conclude that domain walls act as 
effective pinning centers, and pinning becomes stronger 
for a finer and more interwoven domain structure and in 
the presence of static magnetism. 



DISCUSSION 

As can be seen from Fig. [2j the critical current density 
as a function of doping has clear asymmetry and peaks 
in the underdoped regime. One may suggest that the 
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FIG. 2: Top panel. Doping phase diagram of the structural, 
T 3 , the magnetic, T m , transformations as well as supercon- 
ducting transition, T c , in Ba(Fei_ a; Co : z : )2As2 single crystals 
0, [2(J. The T s (x) and T m (x) split with the increasing x. Red 
line inside the superconducting dome shows 0.7T C at which 
the critical current (bottom panel) was evaluated. Inset in 
the top panel shows doping dependence of the degree of the 
orthorhombic distortion S(x), taken at 7 K, base temperature 
of X-ray measurements. Bottom panel shows the evolution 
of the superconducting critical current density j Ct max (stars). 
The critical current was evaluated at a fixed reduced temper- 
ature T/T c (x)=0.7 (a dotted line in the top panel) and in the 
field corresponding to the maximum of the fishtail magneti- 
zation, see Fig. [3] 



observed enhancement is simply due to more robust su- 
perconductivity at the optimal doping. To refute this 
argument, let us consider two samples with similar T c on 
two sides of the dome. For example, for x — 0.054, T c = 
19.2 K and j c = 1.36 x 10 5 A/cm 2 , whereas in overdoped 
x = 0.086 with even higher T c = 20.5 K, the critical 
current is substantially lower, j c = 4.8 x 10 4 A/cm 2 . 
Thus another reasons for the enhancement of j c should 
be looked for. 

Comparison of top and bottom panels of Fig. [2] sug- 
gests that the critical current peaks in the range where 
structural domains (1) become finely spaced, (2) their 
density increases, (3) they form interwoven patterns; (4) 
domains are more likely to be associated with both struc- 
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FIG. 3: Left panel. Evolution of the magnetization M vs. 
H loops with doping at T/T c = 0.7. The magnetization was 
converted into current density using Bean model [lj] to allow 
for a direct comparison of samples with different sizes. The 
data for different doping levels were taken at the same reduced 
temperature T/T c (x)— 0.7. Arrow shows the definition of H m 
(in this case for x =0.061) used to plot j c ,max in Fig. [2] Right 
panel. The temperature dependence of the critical current 
density for sample with a;=0.058, line shows a power-law fit. 



tural and magnetic order. As we will argue, ALL these 
features are favorable for pinning. First, the pinning 
strength is maximal when the size of the pinning centers 
(domain walls and their intersections) becomes compa- 
rable to the size of vortex cores. Second, volume pin- 
ning force increases with the increase of the density of 
the pinning centers. Domain pattern becomes finer with 
x presumably due to lower energy cost of the domain 
wall formation with decreasing distortion 5. Therefore, 
by moving towards finer and denser domain structure 
we both increase the pinning potential (more interwo- 
ven structures) and introduce more pinning centers per 
unit volume. Domains in Ba(Fe i _ a; Co a ;)2As2 propagate 
throughout the entire sample and, for x close to opti- 
mal, form a dense maze of intersecting domain bound- 
aries. A more frequent occurrence of interwoven patterns 
is also caused by the decreasing 5. The areas of cross- 
ings are characterized by strong lattice deformation and 
their formation is energetically unfavorable as compared 
to parallel domain boundaries. Such natural linear de- 
fects extend along the c-axis, similar to heavy ion tracks. 
The nature is generous to at least 122 family of pnictide 
superconductors to create the defects of this type close 
to the optimal doping conditions where T c is the highest 
and the superconductivity is most robust. The contribu- 
tion of the magnetic order may also be very important for 
pinning, as suggested by the data in the bottom panel of 
Fig. [2 Because magnetic order in the planes requires re- 
duced symmetry, orthorhombic distortion inevitably in- 
troduces disruptions of the long range magnetic order. 
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Magnetic modulation itself creates strong modulation of 
the superconducting order parameter, acting as pinning 
centers. This modulation can be enhanced by the magne- 
tostriction, as was shown by decoration imaging of vortex 
lattice in magnetic RN12B2C 0, 0]. 

The critical current density in Ba(Fei_ a; Co a ;)2As2 
reaches maximum at a composition with x=0.058. In 
addition, the critical current has a strong temperature 
dependence, as shown in the right panel of Fig. [3] An 
extrapolation to liquid helium temperatures gives values 
above the "technologically useful limit" of 1 MA/cm 2 . 
This is almost 3 times higher than in the optimally doped 
crystals with cc=0.074 [l5|,[ll|. Considering that the sam- 
ples of Ba(Fei_ x Co :c )2As2 in this range have upper criti- 
cal fields as high as 60 T [a. Hoi [Tll| . critical temperatures 
of above 20 K, low anisotropy 11, [I2> [H| and, as 
shown here, high intrinsic critical current densities, these 
materials can be considered as good candidates for ap- 
plications. 
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